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S : ABSTRACT 

In the minimal supersymmetric model (MSSM), the upper limit of the lightest Higgs boson 
mass, m/jo, depends strongly on the top quark mass, irtt. We have computed the dominant two- 
loop radiative corrections to this upper limit of m^o in order to eliminate large uncertainties 
due to QCD and vn^ corrections. It is shown that the QCD corrections significantly reduce the 
one-loop corrections. As a result, the SUSY parameter space accessible to LEP experiments is 
significantly increased. 



In recent years supersymmetric theories have become maybe the most popular alternatives 
to the standard model (SM) of elementary particle physics. In the minimal supersymmetric 
extension of the SM (MSSM) the Higgs sector contains only the two doublets, Hi and H2, 

required to give masses to up and down type fermions with all the quartic couplings related to 

2 
the gauge couplings. This leads to various restrictions among the Higgs masses and couplings . 

The most important consequence is the existence of a well defined tree-level upper limit for the 

mass of the lightest Higgs boson 



m^o < m I cos 2/3 1 < rriz , with m = min{?T2z, rrij^o} . (1) 



However, it has been shown recently that radiative corrections can significantly alter this pre- 
diction . In particular, the experiments at LEP200 may not be able to detect h'^ or rule out the 
MSSM. The region in the SUSY parameter space, than can be ruled out at LEP experiments 
depends crucially on the top quark mass, mt. However, there is a significant uncertainty in mt 
due to large QCD corrections {e.g. the running mass differs from the pole mass by (9(7%) ). 
Eliminating this uncertainty by defining mt as the pole mass at the one-loop level requires an 
explicit two-loop calculation of m/jo . 

The case tan j3 -^ 00 is particularly interesting because here the tree-level constraint rrif^o < 

rUz is saturated as long as m^o ^ 171^ and thus we expect this case to yield the maximum Higgs 

mass. Note that m^o has another maximum in the limit tan/5 = {i.e., V2 = 0). However, the 

constraint that the Yukawa couplings do not develop a Landau-pole at high energies, together 
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with the experimental lower bound on mt > 113 GeV (95% CL) require that tan/5 ^ 0.5 . 
On the other hand, f3 ~ tt/2 is theoretically very favorable since it would explain the large ratio 
of mt/mi, {e.g. grand unified theories based on SO(IO) predict tan/5 = mt/mi,). 

The numerical analysis of the one-loop corrections ' ' shows that the dominant contribu- 
tions to m^o come from the top-stop sector due to an gfml dependence {gt is the top Yukawa 
coupling). Thus we expect the dominant two-loop corrections to be the contributions propor- 
tional to gfrnif and g'^g'^m.'^ {g^ is the QCD gauge coupling). These terms can be obtained most 



easily in the approximation where we set g and g' to zero. In this case the tree-level potential 
of the lightest Higgs doublet reduces simply to 

V^W = mlHlH2 . (2) 

(Note that here m2 = 771^0 = at tree-level.) After shifting the neutral CP-even component by 
the vacuum expectation value (VEV) [if 2 -^ (^^ + vo)/\^] we obtain 

V!^^ = th' + ^,ml,ih'f + Oiihy] , (3) 

where toj^oq = m^ and t = vqiti^. If we eliminate m2 in favor of t we end up with the relation 
^/i^O = t/vQ. Remember that m/jOQ is still an unrenormalized parameter which has to be 
expressed in terms of a physical observable. This is the physical mass of /i^ (without a subscript 
0) which is identified in the usual way as the pole of the propagator 

^lo = ^hOQ + Re AhOho{mlo) . (4) 

The self-energies and one-point functions (tadpoles) are defined in fig. 1. Note that hP is a 
stable particle at tree-level (the decay of h into gluons is induced at the one-loop level, and 
we can assume that m^o < 2rmn{mt,m^^^mi^}). Thus, eq. (4) is valid even at the two-loop 
level. We now demand that vq is the true VEV to all orders in perturbation theory. This means 
that the tadpoles corresponding to a Higgs field disappearing into the vacuum are absent. This 
reads t + A/jo(0) = or rn^QQ = m^ = —Af^o{0)/vQ. Thus the two-loop renormalized Higgs mass 
can be written as 

m^o = AhOho{m^^o) . (5) 

VQ 

The calculation of the two-loop diagrams can be simplified considerably by using the fact 
that scalar one-particle irreducible n-point Green functions with zero external momenta can 
be obtained as derivatives from the effective potential, Vgfj, with respect to the corresponding 



fields, e.g. 



mlvo + Aho{0) and „ " = mg + A/,ofeo(0) . (6) 



dt>o d^fO 

Furthermore, it is easy to show that 

(— - ^-] Kff = 4^ (^ ' Veff , (7) 

where Vgfj on the right hand side is a function of the masses and couphng constants, but not 
of vq. Thus we arrive at 
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m^o=4^ — ^ V^s + AhOho{mlo) - AhOho{0) . (8) 

Vq \ am^Q J 

The advantage of using the effective potential is that we only have to compute the difference 
A/jO/j()(m^o) — A/jo/jo(0). Note that in a two-loop calculation it is sufficient to compute only the 
one-loop contributions to the Higgs self energies since m/jo here is generated at one-loop. In 
particular, we have to keep the 0[e) terms of the one- loop radiatively corrected Higgs mass, 
m/jo, on the right hand side of eq. (8), since it occurs in a divergent expression. In addition, 
we have to compute the effective potential to second order [\4fj = V^^' + V^^' + V'^-*]. The 
results of this calculation are presented in Appendix B. In the terms derived from V*-^-* we are 
again allowed to replace the unrenormalized quantities by the physical ones. However, V^^' 
depends on the unrenormalized top quark mass, rriiO) the left- and right-handed top squark 
masses, M^ q {P = L, R; we ignore the possibility of L-R mixing due to trilinear Higgs-squark- 
squark coupling. At, which yields only non-logarithmic and thus negligible corrections), and 
the VEV, Vq, at the one-loop level. Thus, we have to express the bare quantities in terms of 
physical quantities to first order. The unrenormalized masses can easily be expressed in terms 
of their pole masses. For the renormalization of the VEV we choose the tree-level relation 
VQ = 2-^/^G^^^^, where G^ = (1.166 39 ± 0.000 02) x 10"^ GeV^^^ jg ^j^g fermi-constant (note 
that the only contributions to the /i-decay come through the W self-energy). Therefore, our 



renormalization conditions are 



vl^2 


-1/2^-1 


(1- 


^m = 


m2(l- 


a), 


<o = 


= <(! 


— a 



(9) 



P = L,R, 



where we have introduced the abbreviations 



a 



^" Up ' ""- m2, ' (10) 

a ^ Re [S^(«) + mu) - Sf («) - Sf (..)] , 

and we have defined u = m^ and Up = M? (P = L, R). It is now straightforward to derive an 
expression for the two-loop radiatively corrected Higgs mass in terms of physical quantities 

("^ft")2LP = ^ (1 + Ow - 2a) 

+ AV2Ncu'^Gf,K^^ I (2a - a^ - a^) 

+ e [(a - ^aw - |a/i) ((In-u^)^ + (In-u^)^ - 2(ln-u)^) 

-l-a^ln-Uj^ + a^ln-u^ — 2aln-u] > 

+ aV2u^G^ (^^ \/(2) ^ AhOho{h) - AhOho{0) , 
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where k = IGvr^ and we have defined 

dA/,Ofto(p2) 



ah = 



dp^ 



(12) 



p'^=h 



We have denoted the one-loop radiatively corrected squared Higgs mass in the approximation 
introduced above by 

h = sV2NcU^G^K-ho . (13) 

Here we have defined to = ln{M^^i^Y / m^) where we found it convenient to parameterize our 
numerical results by Mg^gy — ^i ^t ^"^^ ^ = ^t l^t ■ '^^^ analytic results for the relevant 
self-energies and the expression for the one-loop and the two-loop effective potential are given 



in Appendix A and B, respectively. Here we will simply present an approximate result in the 

case of heavy, mass-degenerate superpartners M^ = M? = M~ ^ m^ 

(^IOapp = h{l- QqICfk-^ [to + 2] 

(14) 
+\g'iK'^ [9to + 18 + 7iVe + (15 - tt^) t^^] ] . 

The values of the higgsino mass parameter, /x, and of the heavy Higgs doublet, m^o, turn out 
to be irrelevant. 

We will now proceed to compare our diagrammatic two-loop result with the second order 
terms obtained by a renormalization group (RG) approach. Here we reintroduce the Higgs 
self coupling. A, in the potential of eq. (2) which we treat as a running low energy effective 
parameter rather than a bare parameter. In this case we find in general that A 7^ at any 
energy scale ^/s < MgusY even in the case g = g' = 0. The RG improved Higgs mass can be 
written as 

where the running quartic Higgs self coupling evaluated at the scale s = mj^o is 

Inrnj 

A= j dt/^A, (16) 

In M|usY 
By using the f3 function to two-loop order, (3x, given in ref. 10 and 11 we obtain the leading and 
next-to-leading log radiative corrections to the Higgs mass summed to all orders in perturbation 
theory. In order to obtain an analytic result we will solve the RG equations (RGE) iteratively. 
The second order contributions to order g^ and g'^gf are 

(mlo) =h{l- toK-^ [6g^,CF - igf] } 

(17) 
("^ft02/3 = ("^fto)i/3 + ^'^ ^ [^9sCf - 5gt] , 

where the subscript 1/3 (2/9) indicates the use of one- loop (two-loop) P functions. It is instructive 
to compare eq. (17) to eq. (14). We see that the coefficients of tg obtained by solving the one- 

Q 

loop RGEs to second order are in agreement as they should be. However, the terms linear in 
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to obtained by solving the two-loop RGEs to first order are different. This is not surprising 
since the RG approach does not include any threshold effects. For example, the QCD threshold 
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corrections to mt are 

a. 



run / 1 , /^ s 

mt = ml ' 1 I ^^ 

n 



(l + Cf- 

where the running top quark mass is tti™'^ = 2^^'^G^ fi'i(^^t) and we use Og = 5's/47r = 0.11. 
We now investigate the origin of the m.f dependence of the one-loop corrections to m/jo more 
closely. From eq. (15), (16) and (18) we obtain 

mlo oc gfG-ho oc mfG^to (l - 4Cf ^) . (19) 

The dg-term in eq. (19) coming from finite threshold corrections to m,t are of the same order 
as the next-to-leading log term of eq. (17). If we add both terms together we indeed recover 
the coefficient proportional to g'^ in eq. (14). 

We will now present the numerical results of our two-loop calculation. It is contrasted 
with the second order leading log result and next-to-leading log result obtained from a simple 
analysis of the one-loop RGEs and two-loop RGEs not including QCD threshold corrections 
[eq. (17)]. We set /x = m^o = 200 GeV, r = 1 and M? = M? = M^^^y - m^. In fig. 2 we 
present the two-loop contribution to the Higgs mass obtained by different methods. Presented 
are the results from our diagrammatic two-loop calculation (solid curve), the second order 
contributions from a RG approach using one-loop and two-loop (3 functions (dotted curve and 
dot-dashed curve, respectively) and result of our approximate expansion (dashed curve). We 
have parameterized the results by a linear mass shift defined as 

(Amj^o)^ = ^^^ , X = 2LP, ip, 2P, APP . (20) 

The contributions to order gfm."^ (Yukawa) and gtgl'rn't (QCD) are plotted separately. The 
dependence of Am/jO on the remaining parameters, r, M^, /i and mj/o changes the result by 
less than 1 GeV as long as Mg, ji < 3Msusy and 0.1 < r < 10. 



Sofar we have focused our attention on the case of large tan (3 where one obtains the upper 
hmit for m^o. We have seen that the diagrammatic result differs significantly from the RGE 
results [eq. (17)]. However, the agreement with our large MgusY expansion [eq. (14)] is excellent. 
Furthermore, the two-loop result is dominated by the QCD corrections which are even larger 
than those one-loop electroweak corrections which are not enhanced by powers of mt and thus 
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never exceed a few GeV . 

In recent years it has become standard in experimental analyses to only include the domi- 
nant mf term at the leading log level via the CP-even mass matrix in electroweak eigenstates 

„ 9 / sin /5 — sin /5 cos (3 

\ — smpcosp cos p 



(21) 

/ nriG^ f-i — Gin f-i nntG f-i \ 



cos^/9 — sin/3cos/3 



* - sin /9 cos (3 sin f3 + em. 



z 



where e = h/ sin /3. Note that at this level our two-loop calculation can be directly generalized 
to the case of arbitrary tan (3 by using 






(22) 



1 + ^ (to + 2) 
n 



27r2sin2/3^ 
where to = ln{Mf Mf /rri^). In fig. 3 we compare this approximation (solid curve) with the 



one- loop result, ^/■m'^ + h, (dotted curve) and the two-loop result, x/m^ + (^^o)2LP! (dashed 
curve) in the limit of large tan/?. We plot rrif^o as a function of MgusY for ^^t = 150 GeV 
and r = 1, and we set mjjo = fi^ = M? = M? = Mg^sY " ""^i- ^^ ^^^ ^^^^ ^^ ^^^ ^^^ 
TeV region the two-loop corrections are comparable to the one-loop result in magnitude due 
to large logarithms. Those terms have to be summed over using the RGBs. We have moved 
the QCD corrections into the denominator which is in better agreement with the RG improved 
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result . In particular, the approximation of eq. (21) and (22) increases monotonically even for 
very large tg. 

In fig. 4 we show contours of constant m/jo = 60, 80, 100, 110, 120, 130 GeV in the tan/?- 
MgUSY plane in the approximation from eq. (21). We compare the results including QCD 



corrections (solid curves) and without QCD corrections (dashed curves) in the hmit mj^o ^ rn^ 
where we obtain the upper hmit of mf^o. 

We have shown that the QCD corrections to m^o are large and negative and thus lead to 
a reduction of the one-loop result by (9(30%) for MgusY = 1 TeV. Consequently the region in 
the SUSY parameter space accessible at LEP experiments increases significantly. The effects 
are particularly important in the case tan/3 = 1 where the tree- level Higgs mass vanishes. The 
implications for the large tan/? regime depend on whether it is possible to detect a Higgs boson 
with a mass, -m/jo < rriz- Consider e.g. the scenario where mt = 150 GeV and the experimental 
lower limit of m^o > 110 GeV. Then the QCD corrections shift the lower limit of MgusY form 
670 GeV to 1.3 TeV . 

Acknowledgements: We would like to thank H. Kiihn for suggesting this project. 



APPENDIX A 

In this appendix we present the results of the one-loop self energies of top quark, top 
squarks, the Higgs boson and W boson in the approximation that g = g' = Ajj = 0. We 
have done all the calculations in dimensional reduction in order to preserve SUSY (i.e., we 
set the number of space-time dimension to n = 4 — 2e, but keep the number of components 
of all other tensors fixed) . We have suppressed the arbitrary renormalization scale. The 
results are presented in terms of the following integrals formulated in Minkowski space using 
the Bjorken-Drell metric 



w .A . . d^A; 1 

AQ[m ) = —IK 



{2tiY k^ -m'^ + iS' 

Bo{q'^,ml,ml) = -in / 777^77^ 2 , -.wd , \2 2 , -xi ' 

J [2t[)^ [k^ — m| + io)[[k + qY — mg + io\ 

q^Bi(q^,ml,ml) = -m / 77^-^772 2 , -nw ^1 / ^2 2^~^ ' (^'^^ 



The top quark self energies are 



AcE^ (u) = —g^ Bi{u, h, u^) + Bi{u, u, 0) 

+ g^CF[-2Bi{u,g,u^) 

+ {1 + OBo{u,0,u) + 2Bi{u,0,u)] , 

i^^t (u) = -gt Bi{u,h,Uj^) + Bi{u,h,di^) 
+Bi{u,u,0) + Bi{u,0,0)] 
+ 9sCF[-'2Bi{u,g,Uji) 
+ {1 + OBo{u,0,u) + 2Bi{u,0,u)] , 

Kmu) = g!CF{3 + OBo{u,0,u). 



(A.2) 



(A.3) 



(A.4) 
(A.5) 



The top squark self energies are 



+ [4uiBo{u^, 0, Uf) - Ao{uf) 

+2 [{u + g- u^) Bo{u^, g, u) + Ao{g) + Ao{u)]} 



9t i [Ao{u) + Ao{h) + {u + h- U];)Bq{uj^, h, u) 
-2uBo{uf, Uf, 0) - uBq{uf, di, 0) - hBQ{u^, u^, H) 

-Aq{Uj^)} , 



(A.6) 



+2 [{u + g- Uf) Bo{uf, g, u) + Ao{g) + Ao{u)]} 



+ 9t 



Aq{u) + AQ{h) + {u + h- Uf)Bo{uf, h, u) 
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+ y^oCh) + Ch- u^)Bo{uj^, h, 0) 

-2uBo{u^, M^, 0) - hBo{uj^, dj^, H) - hBo{uji, u^^, H) 

-Ao(«J-Ao(4)|. (A.7) 



The Higgs boson self energy is 

KAhn,o{h) = Ncg^ {[(4m - h) Bo{h, u, u) + 2Ao{u)] 

-2uBo{h, uj;^, Uj^) - 2m5o(/i, m^, Uji) (A. 8) 

-M^l)-Mur)} . 
The W self energy is 

nA^wiO) = N^mlg^ f 1 + i - lE« + ^^^ - ^ In §^ V (A.9) 

where we abbreviate the gluino mass, the left-handed bottom squark, the Higgsino mass and 

f? h = jj,'^ and H = m'jjo 



the mass of the second Higgs doublet mass hy g = M? dj^ = M? , h = fi'^ and H = m'jjo, 



respectively. 



APPENDIX B 



Here we present the calculation of the Feynman diagrams contributing to the effective 
potential. We have done all the calculations in dimensional reduction to preserve the super- 
symmetric Ward-identities. The QCD gauge sector has been calculated in arbitrary R^ gauge. 
The results for the one-loop and two-loop effective potential are presented in terms of the 
following integrals 



K{m' 






and 



K J{mi,ni2) = AQ{mi)AQ{m2) 
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It I I I l\ 



(B.2) 



Furthermore, we find it convenient to introduce 



L{a, b, c) = J{c, h) — J{a, b) — J{a, c) — {a — b — c)I{a, b, c) . 



The one-loop effective potential is 



V-d) = AT, [2Kiu) - K{ul) - K{u^)] 



The relevant QCD contributions to the two-loop effective potential are 



v; 



(2) _ „2/Ar2 



gt{N', - 1) {[-2uI{u,u,Q) + J{u,u)] 



+ [L{u^,g,u) + L{u^,g,u)] 
and the second order top Yukawa contributions are 



K 



(2) 



Ncgt \ L{0, u, u) + L(0, 0, u) + L(m^, h, u) 



+L{u^, h, u) + L{di^, h, u) + L{u^, h, 0) 

- ul{u^, u^, 0) + ul{d^, u^, 0) + m/(m^, Ur, 0) 



+ 



J{ul,Uji) + Jid^^up.) 



(B.3) 



(B.4) 



(B.5) 



(B.6) 



The derivative of V^'^' = \4 + V^ in eq. (11) is obtained by using the chain rule and 

duL duR Ml dg 

du du du du 

We have computed the second derivative of V^'^' numerically and analytically. Furthermore, 
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we have checked both analytically and numerically that the right hand side of eq. (11) is 

independent of e and the renormalization scale. 

Remark : We have done the two-loop calculation in dimensional reduction. Operationally this 
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was done using dimensional regularization plus the contributions from the so-called e-scalars . 
The latter contribution is divergent due to the 1/e^ pole from the two-loop integral. Without 
it eq. (11) is not finite. 
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FIGURE CAPTIONS 

1) The definition of tlie gauge boson top quark and scalar self energies, and the Higgs 
tadpole. The left and right handed projection operators are R, L = {1 ± 75)/2. 

2) The two-loop shift in the Higgs mass, (^Amj^o), as a function of the top mass. We present 
the terms to order gfm^ (Yukawa) and g^glml (QCD) separately. 

3) The radiatively corrected Higgs mass, m/jo, as a function of MgusY- Included are terms 
to order gfml and g^glml. 

4) Contours of constant Higgs mass, m^fi, in the tan/3-MsusY plane. We take mt = 150 GeV 
and consider the limit m^o ^ ^^z- 
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